Diamagnetic Phase Transition and Phase Diagrams 

in Beryllium 
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The model of diamagnetic phase transition in beryllium which takes into account the quasi 2- 
dimensional shape of the Fermi surface of beryllium is proposed. It explains correctly the recent 
experimental data on observation of non-homogeneous phase in beryllium at the conditions of strong 
dHvA effect when the strong correlation of electron gas results in instability of homogeneous phase 
and formation of Condon domain structure. 

PACS numbers: 75.20.En, 75.60.Ch, 71.10.Ca, 71.70.Di, 71.25.-s: 71.25. He; 75.40-s; 75.40.Cx. 



I. INTRODUCTION 



It is well-known |lj that the instability in electron gas 
due to magnetic interaction between conduction electrons 
in diamagnetic metals at high magnetic field and low 
temperature leads to formation of non-uniform diamag- 
netic phase, so-called Condon domains, which is usually 
realized as a stripe-domain structure for plate-like sam- 
ples. The diamagnetic domains were observed in beryl- 
lium by magnetization measurements [3 an d m silver 
by nuclear magnetic resonance (NMR) [3j and by Hall 
probes spectroscopy 0. They have been observed in 
beryllium, white tin, aluminium and lead by muon spin- 
rotation (/xSR) spectroscopy 0. The above-mentioned 
instability of an electron gas is called diamagnetic phase 
transition and is extensively studied, both theoretically 
and experimentally [t|-[l2|. 

Recent progress in experiments on observation of Con- 
don domain structure in silver and beryllium 0-0 
provides a natural stimulus towards a more detailed un- 
derstanding of the properties of strongly correlated elec- 
tron gas at the conditions of dHvA effect. Undoubt- 
edly, silver and beryllium are the most popular met- 
als in experimental investigation of electron instability 
leading to formation of Condon domain structures. The 
direct observation of Condon domains in silver by Hall 
probe spectroscopy shows that the temperature and 
magnetic field dependencies of non-uniform phase are in 
good agreement with the theoretically predicted phase 
diagrams (see, e.g. estimated on the basis of the 

Lifschitz-Kosevich-Shoenberg formalism in the first har- 
monic approximation . However, in case of beryllium, 
the standard approach, which gives the satisfactory re- 
sults for free or almost free electron gas, fails to describe 
the diamagnetic phase transition, and, leads, in particu- 
lar, to the underestimated values for critical temperature 
T c and magnetic field H c , which is inconsistent with the 
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experimental data 0-0- Thus, for proper calculation of 
amplitude of the dHvA oscillations and constructing the 
phase diagrams, the correct topology of Fermi surface has 
to be taken into account. Recent experimental observa- 
tion of irreversible effects in beryllium by Hall probes in 
dc field and standard ac method with various modulation 
levels, low frequencies and magnetic field ramp rates 0, 
as well as application of non- linear techniques 0, which 
reveals giant parametric amplification of non-linear re- 
sponse in Condon domain phase, offer the possibility to 
construct the diamagnetic phase diagrams for beryllium 
and compare it with theoretical predictions. 

The Fermi surface of beryllium was investigated very 
carefully in the past . Paradoxically, beryllium has the 
simplest Fermi surface, because it differs essentially from 
the free electron model. The Fermi surface of beryllium 
consists only of the first and second zone monster ('coro- 
net') and the third zone 'cigar'. It is well-established 
that the dHvA oscillations originates from the three max- 
imum cross-sections of 'cigar' ('waist' and 'hips') which 
are characterized by a very small curvature. The small 
curvature of the cylinder like Fermi surface explains rel- 
atively high amplitude of dHvA oscillations in beryllium 
at magnetic field H applied parallel to the hexagonal 
axis, e.g. H || [0001], in comparing to silver jj], where 
the small amplitude of dHvA oscillations is explained in 
the framework of spherical Fermi surface [l| providing 
the reasonable agreement between estimated phase di- 
agrams | l2j and experimental data 0. A new experi- 
mental results 0-0 on observation of Condon domain 
phase in beryllium stimulate the further development of 
the theory. Undoubtedly, for correct explanation of ex- 
perimental data |9-|3 the anomalously low curvature of 
the 'cigar' like part of Fermi surface of beryllium near 
the extreme cross-sections ('waist' and 'hips') has to be 
taken into account. The modeling of cigar-like part of 
the Fermi surface of beryllium |13| by a cylinder, similar 
to 2D electron gas, describes correctly the diamagnetic 
phase diagrams at low range of quantizing magnetic field 
(< 3T ), but results in the essentially overestimated val- 
ues of critical parameters (temperature T c and magnetic 
field H c ) at higher values of applied magnetic field. In 
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particular, the models based on approximation of rele- 
vant Fermi surface sheets by cylinder predict the exis- 
tence of inhomogeneous phase at the values of external 
magnetic field H ~ 10 — 12T contradicting to the ex- 
perimental data 0-0 that show the disappearance of 
Condon domain structure above ~ 6T at the values of 
Dingle temperature Tjj w 2.2K. Moreover, the approxi- 
mation of Fermi surface of beryllium by cylinder 13] fails 
to explain the observed beatings in the dHvA oscillations 
0i0|"0> which are the result of the interference of the 
signals from two different extreme cross-sections of Fermi 
surface with close fundamental frequencies. Solt 
shows the serious disagreement of /iSR data in beryllium 
with the phase diagrams, calculated in the framework of 
standard Lifshitz-Kosevich formula, and necessity to take 
into consideration the actual 3D Fermi surface geometry 
for beryllium, first, the existence of two different cross- 
sections ('waist' and 'hips'), and second, the low curva- 
ture of Fermi surface near extreme cross-sections. The 
model developed in 0-0 is more realistic; it describes 
correctly the beating effect in dHvA oscillations in beryl- 
lium, and it is consistent with the previous experiments 
on observation of Condon instability at low range of 
quantizing magnetic field (till ~ 3T). Unfortunately, 
the high values of critical magnetic field ~ 10T for non- 
homogeneous phase, estimated in the framework of model 
, contradict to recent data on independent observation 
of Condon instability in beryllium 0-0, which demon- 
strate the existence of non-uniform phase in essentially 
narrow interval of quantizing magnetic field (r* — 6T) 
depending on estimated Dingle temperature Tr> ~ 2.2K. 
In present paper we develop the model of slightly corru- 
gated cigar-like Fermi surface for beryllium and calculate 
the diamagnetic phase diagrams T-H at different Dingle 
temperatures Tjj, The estimated phase diagrams are in 
good agreement with the available experimental data on 
observation of diamagnetic instability in beryllium 0-0. 

The paper is organized as follows. Section I is an 
introduction. In Section II we consider the model of 
slightly corrugated cylinder Fermi surface of beryllium. 
The model is based on reliable knowledge of the dimen- 
sions of the relative Fermi surface sheets 0, as well as 
its qualitative nature. In Section III we calculate the 
diamagnetic phase diagrams for beryllium and compare 
the theoretical results with recent observation of electron 
instability in beryllium 0-0. Section IV is conclusions. 
Section V contains acknowledgements. 



II. MODEL 



The oscillatory part of the free energy density which 
is responsible for dHvA effect can be written as follows 

0-0 



where 
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(e/xpff)- 
(27r 2 c) 2 m 



-R(T,^H,T D ) [ dk z cob[27tX(A,)], (1) 

c J 



R(T,^H,T D ) 



A( Mo ff)Texp [-X(fioH)T D ] 
sinh [X(^ Q H)T] 



(2) 



Here, e is the absolute value of electron charge, c is 
the light velocity, ks is the Boltzmann constant, m c is 
the cyclotron mass, X(fiaH) = 2ir 2 kBm c c/ 'eh^H and 
Td = fr/2TtkBT is Dingle temperature which is inversely 
proportional to the scattering lifetime t of the conduc- 
tion electron. The cross-sectional area of the Fermi sur- 
face A(fi,k z ) at k z in k-space is related to the quantity 
X(k z ) by the relationship 



X(k z 



chA(\ik z ) 
2irefioH 



7, 



(3) 



where fj, is the chemical potential and 7 is the phase cor- 
rection (typically, it equals to 1/2). 

The integral in Eq. is a Fresnel-type. Its major con- 
tributions come from regions where the phase is station- 
ary 0, e.g. the cross-sectional area A(fi,k z ) has maxi- 
mums or minimums. The usual procedure of calculation 
of such an integral consists in expanding of electron orbit 
area A(fi,k z ) about the extreme points and summation 
of the contributions from different extreme cross-sections 
of Fermi surface 0. In case of beryllium the standard 
procedure requires considerable modification due to van- 
ishing low curvature of extreme cross-sections of 'cigar', 
when the higher order terms in expansion need to be 
considered. The additional disadvantage of standard ex- 
pansion of A(n, k z ) about the extreme points arises from 
the fact that neither curvature factor nor its derivatives 
are available from the experiment. The last circumstance 
implies the necessity of development of model represen- 
tations of the Fermi surface of beryllium. 

Our consideration is based on the following model rep- 
resentation of electron orbit area for the 'cigar' like sheet 
of Fermi surface of beryllium relevant for observed dHvA 
oscillations: 



A(k) 



1 + 6, 0<|fc|<l-/3 

1 + Ssma(k - 0.5), (3 <\ k \< 1 - f3 

1 + 5, 1-/3<| fe |< 1+/3 

l-(5sina(fc-0.5), 1 + (3 <| k |< 2 + (3 



(4) 

where k = k z j A is reduced wave vector, a — 7f/(l — 2(3), 
(3 = k c /A and A = 0.217 1/A° is the distance in recip- 
rocal space between two extreme cross sections of 'cigar'. 
The quantity S = (Ah — A w )/2Ao characterizes the dis- 
crepancy between maximal (Ah) and minimal (A w ) cross 
sections of Fermi surface, and Ao — (Ah — Aw)/ 2 is aver- 
age cross section area. According to Eq. (0} the extreme 
cross sections of 'cigar' are characterized by the same cur- 
vature which is consistent with the data on Fermi sur- 
face of beryllium 0. In calculation of phase diagrams 
one must take into consideration the transition region 
between 'waist' and 'hips'. We use the approximation 
of the transition regions by trigonometric functions (see, 
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Eq. i|Hl). defined in corresponding intervals through the 
adjustable parameter < (3 < 1/2. Although, the Fermi 
surface of beryllium was investigated very carefully in 
the past 0-0 > still there is lack of data concerning the 
transition between extreme cross sections. We show that 
with suitable choice of this parameter the satisfactory 
agreement between calculated phase diagrams and ex- 
perimental data can be achieved. Before calculation of 
the phase diagrams we emphasize that our choice of the 
representation of the Fermi surface sheet of beryllium in 
the form of Eq. (0} is governed by the experimental facts 
0-0 which indicate on the negligible and equal (or al- 
most equal in the accuracy of experiments) curvatures of 
extreme cross sections. The last circumstance, e.g. negli- 
gibility and equality of the 'waist' and 'hips' curvatures, 
simplifies our task and allows us to neglect the curvatures 
in the vicinity ~ 2/3 of the extreme cross sections with 
proper choice of the functions describing the transition 
between two extreme cross sections. 



III. RESULTS 

The oscillating part of the diamagnetic susceptibility 
X — —dM = —d 2 fl/(dB) 2 is obtained by differentiat- 
ing Eq. Q twice with respect to the magnetic induction 
B = hqH + M. Calculating the integral in Eq. © with 
taking into account the representation (@J) and keeping 
the leading terms only (see, e.g. Q), one can arrive to 
the following expression for the susceptibility 



X = a cos (l 2 Ao + ip), 



(5) 



where 



« = t A( f" 1 I Q(^oH) | R(T, noH, T D ), (6) 

is reduced amplitude of dHvA oscillations, e.g. the ratio 
between the amplitude of dHvA oscillations and their 
period. 

The complex function 

Q = Qi + Q 2 =| Q exp jtp, 
Qi =3/3cosx + 2(l-2/3)J (x), Q 2 =P$inx (7) 

where x = 1 2 AqS, depends on applied magnetic field 
through parameter I 2 — ch/e^H. Here, Jo(x) is Bcsscl 
function of the first order. 

When the amplitude of dHvA oscillations become high 
enough (a > 1 Eq. ©) the magnetic interaction between 
electrons results in thermodynamic instability of uniform 
phase and diamagnetic phase transition into non-uniform 
phase with formation of Condon domain structure oc- 
curs. The condition 



a(T,[i H,T D ) = 1 



(8) 



defines the critical curves on the plane T — H at differ- 
ent Dingle temperatures Td- The Dingle temperature 
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FIG. 1: Shown is the field dependence of the Dingle tem- 
perature Td = Td(hoH) at T — 0. Depending on the impu- 
rity of the sample, there is a set of intervals in the applied 
magnetic field with possible existence of the Condon domain 
structure. The dash line corresponds to the value of To = 2K, 
it crosses the curve To = Td(^oH) in discrete number of 
points which divide the whole range of magnetic field into al- 
ternating intervals. With increase of the field H the uniform 
phases are replaced by non-uniform ones with the periodicity 
P(l/fi H) = 1/AF. In calculation we put (3 = 0.15. 



Td characterizing the quality of the material can play a 
crucial role in observation of electron instability. Thus, 
for comparing the theoretical results with experiment it 
is important to know the values of Td , for which the 
Condon domain structure can be observed in principal 
at given value of magnetic field H . It is convenient to 
consider the function Td = Td(^qH), defined by the 
Eq. © at T = OK. Solving the Eq. © at T = OK, one 
can obtain 



T D (noH) = A(Mo#)m 



(2A h) 2 k 1 
7T 3 ra c (/i #) 



T Q(poH) (9) 



The results of numerical calculations of the function 
T D = Td(hoH) © are shown in Fig.Q 

The function Td = Td{^qH) is characterized by os- 
cillatory dependence on the magnetic field H, which is 
a consequence of the beats resulting from two close fun- 
damental frequencies. It is interesting to compare this 
dependence with analogous one for silver. In case of sil- 
ver [l5j |. at given impurity of the sample, e.g. fixed Din- 
gle temperature Td , the values of external magnetic field 
with possible existence of Condon domain structure be- 
long to the interval H < (T D < H < H > {T D ))- For beryl- 
lium the whole interval of quantizing magnetic field con- 
sists of the set of the intervals H { <\t d < H < H^{T D )) 
with alternating uniform and non-uniform phases. The 
number of the intervals, where the non-uniform phase ex- 
ists, decreases with increase of Dingle temperature Td), 
collapsing at some value T D max ^ , depending on parameter 
(i. Thus, T^""" w 3.1A' for f3 = 0.15. For T D = 2K the 
width of the interval of the existence of the non-uniform 
phase, defined as AH 1 = i2> — if < , increases with in- 



creasing the external magnetic field H. One can show 
that the maximums of the function appear periodically 
on the scale of inverse magnetic field with the period in- 
versely proportional to the discrepancy of the two funda- 
mental frequencies, corresponding to two extreme cross 
sections of 'cigar': 



where AF = F h - F w . With F h = 970.9T and F w = 
942.2T we obtain AF = 28. IT. This can be verified 
by experiment on observation of non-uniform phase. 

To construct the phase diagrams for beryllium we put 
a = 1 in Eq. |JBJ. The T — H curves in Fig. Inform the ge- 
ometric place of points which correspond to diamagnetic 
phase transition at different Dingle temperatures. Out- 
side the shell- like surface To = To {T, H) the uniform 
phase takes a place, while the ordered phase is located 
below this surface. At given Dingle temperature To the 
field dependence of the critical temperature T = T(fioH) 
has maximums. The analysis of the phase diagrams esti- 
mated in the framework of Eq. 0} shows that the posi- 
tion of the maximums depends entirely on the difference 
between fundamental frequencies, corresponding to two 
extreme cross sections of the Fermi surface, while the val- 
ues of these maximums is defined by the impurity of the 
sample, e.g. the values of Dingle temperature To, and 
also by the parameter (3 < 0.5, which is the character- 
istic of the transition range between two extreme cross 
sections of Fermi sheet. The analysis shows that the val- 
ues of the maximums grows with decrease of the Dingle 
temperature To or increase of (3. The results of numer- 
ical calculations of the phase diagrams for different Din- 
gle temperatures To = 0, 1 and 2.1K are illustrated by 
Fig. [31 which shows that the positions of the maximums 
of the curve T — H are not affected by the Dingle temper- 
ature, but the values of the maximums depend strongly 
on the To ■ The increase of To results in decrease of the 
values of the maximums and their collapses. Thus, the 
last maximum, located at for To = OK, disappears when 
T D grows to 2.1K. 

Resent experimental results on observation of irre- 
versible effects in beryllium Q, as well as giant para- 
metric enhancement of non-linear effects , allows us to 
verify the phase diagrams calculated in the framework of 
the assumptions 10} . Undoubtedly, both of the observed 
effects, e.g. irreversibility |6( and strong non-linear dy- 
namics are closely related to the electron instability 
at given temperature in increasing applied magnetic field. 
This instability results in the phase transition from uni- 
form state to non-uniform one with formation of Condon 
domain structure. Both results 0]-0 indicate on the 
existence of the large scale periodicity of the measured 
effects relative to inverse magnetic field (see, Eq. (|10[l ). 
which is additional one to the usual small scale periodic- 
ity of dHvA effect. Undoubtedly, the large scale periodic- 
ity of above-mentioned effects, measured at given temper- 
ature (T = 1.3K @ and T = 1.2K 0) is the result of the 




FIG. 2: Three dimensional phase diagrams in coordinates T— 
H -T D for beryllium. Above the surface T D = T D (T,^ H) 
the uniform phase exists, while below this surface the non- 
uniform phase (Condon domain structure) is formed. 

periodicity in the appearance of the non-uniform state 
and, therefore, can serve as an instrument for verification 
of the theoretical phase diagrams. Comparison between 
the theoretical phase diagrams calculated on the basis of 
the assumption 01 is illustrated by Fig. 01 In Fig.0{a) the 
calculated diamagnetic phase transition temperature is 
plotted as a function of applied magnetic field for Dingle 
temperature To = 2.2K. In calculation of phase diagram 
we put (3 = 1.15. In Fig.0}6) the amplitude of the third 
harmonic (in arbitrary units), arising in the conditions 
of electron instability, as a function of magnetic field is 
depicted. It is evident from the comparison the theoret- 
ical phase diagram and data 0, measured at T = 1.3-K", 
that non-linear response appears at the intervals of ap- 
plied magnetic field, corresponding to the formation of 
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FIG. 3: Envelope of phase diagrams for beryllium T — H 
(color online), estimated at three different Dingle tempera- 
tures: T D = Q,1K and 2.1K (curves 1) (black), 2 (blue) and 
3 (red) correspondingly. The position of the maximums is not 
influenced by Td, while the values of the maximums depend 
strongly on Td, decreasing with increase of Td- 
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FIG. 4: Comparison computed phase diagrams with the data 
0- Typical critical curve T = T(fioH) at given Dingle tem- 
perature To = 2.2K is shown (a). Figure (6) represents the 
experimental data Q on investigation of the irreversible ef- 
fects in beryllium, which are the results of formation of Con- 
don domain structure. The output signal is shown in arbitrary 
units. There is a good agreement between the intervals of ex- 
istence of the non-uniform phase, estimated in the framework 
of and the appearance of the observed in Q non-linear 
response. 

the Condon domain structure. The inverse large scale 
period AF = 28.7TlT7)l gives a good agreement between 
theory and experiment. Unfortunately, the relevant to 
our calculation measurements in [|| were made only at fix 



temperature T = 1.3K. The measurements of the tem- 
perature dependence of the non-linear response (together 
with the field dependencies) would make it possible to re- 
construct the complete phase diagrams T = T{^lqH) for 
beryllium. 



IV. CONCLUSIONS 

We developed the model of slightly corrugated Fermi 
surface of beryllium with taking into account the negli- 
gible small curvature of the extreme cross sections of the 
Fermi surface sheet relevant to observation of dHvA os- 
cillations. The calculated in the framework of our theory 
phase diagrams for beryllium are in a good agreement 
with available experimental data on observation of Con- 
don domain instability. 

The estimated phase diagrams reveal large scale peri- 
odicity on reciprocal magnetic field. The inverse period 
does not depends on the impurity of the sample and is re- 
lated to the discrepancy of the two fundamental frequen- 
cies AF = 28. 7T Eq. IjlOp , corresponding to two extreme 
cross section of 'cigar' like Fermi surface of beryllium, 
e.g. 'waist' and 'hips'. The magnetic field dependencies 
of recent experimental data on observation of non-linear 
effects in beryllium Q-Q have a simple explanation on 
the basis of the estimated diamagnetic phase diagrams. 
We hope that the theoretical results will stimulate fur- 
ther experiments on investigation of electron instability 
in strongly correlated electron systems. In particular, 
the measurements of the temperature dependencies of 
the non-linear effects in the samples of different impurity 
will allow constructing the complete phase diagrams. 
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